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J> . Abstract 

^^ ' The quantum-classical crossover from the Fermi liquid towards the 

Wigner solid is numerically revisited, considering small square lattice 

C*) . models where electrons interact via a Coulomb U/r potential. We review 

a series of exact numerical results obtained in the presence of weak site 
disorder for fully polarized electrons (spinless fermions) and when the 
spin degrees of freedom are included. A novel intermediate regime be- 
tween the Fermi system of weakly interacting localized particles and the 
ti . correlated Wigner solid is obtained. A detailed analysis of the non disor- 

dered case shows that the intermediate ground state is a solid entangled 
with an excited liquid. For electrons in two dimensions, this raises the 

17—* ■ question of the existence of an unnoticed intermediate liquid-solid phase. 

^^ —1/2 

<**] ' Using the Coulomb energy to kinetic energy ratio r s oc U <x n s , we 

Q , discuss certain analogies between the numerical results obtained as a 

O ■ function of U for a few particles and the low temperature behaviors ob- 

tained as a function of the carrier density n s in two dimensional electron 
gases. Notably, the new "exotic state of matter" numerically observed 
p\. ' at low energies in small clusters occurs at the same intermediate ratios 

5-H t r s than the unexpected low temperature metallic behavior characteriz- 

ing dilute electron gases. The finite size effects in the limit of strong 



disorder are eventually studied in the last section, providing two nu- 
merical evidences that the weak coupling Fermi limit is delimited by a 
second order quantum phase transition when one increases U. 
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"The very simplest form of the theory of the energy bands in metals 
gave for many problems such accurate explanations of often very intri- 
cate properties of metals and alloys that it may well appear superfluous 
to consider extensions of the simple form of the theory" . Those words 
written by Wigner [1] in 1938 come again as an objection against the 
need to develop a more rigourous theory, since the Fermi liquid theory 
(FLT) was improved by Landau [2] in the sense of a perturbation theory 
based on renormalized single-particle excitations and adapted to include 
the effects of elastic scattering by the impurities [3] . The need to go out- 
side conventional FLT for explaining the unexpected two dimensional 
metallic phase [4] discovered by Kravchenko and Pudalov is a subject of 
controversy. On one hand, certain characteristic FLT behaviors [5] seem 
to remain in the vicinity of the metal- insulator transition (MIT), sug- 
gesting that an "apparent" metallic behavior could be the consequence 
of "classical" effects (interband scattering [6], temperature dependent 
screening [7], temperature dependent scattering [8, 9] or classical perco- 
lation [10]). On the other hand, the observation of an unexpected MIT 
is first the result (see for instance Refs. [11, 12]) of new possibilities of 
studying controlled many body systems which are closer to the strong 
coupling limit than the previously studied systems in two dimensions. 
This gave us the motivation to numerically revisit the classic problem 
of the crossover from the weak coupling Fermi limit towards the strong 
coupling Wigner limit for electrons in two dimensions. In this chapter, 
we review our main numerical results. The interest of the information 
given from exact diagonalization of small systems can be questionned, 
but it may have the merit to raise questions which may be relevant for 
explaining the behaviors observed around the "two dimensional MIT" . 

1. WEAK AND STRONG COUPLING LIMITS 

A convenient measure of the electron gas density n s is the dimension- 
less parameter r s 

r s = — L- (1) 
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defined as the radius 1/ \fim s of the unit disk divided by the Bohr radius 
a B = h /(me 2 ). The unit disk encloses an area equal to the area per elec- 
tron of the gas. For a real two dimensional electron gas (2DEG) created 
in a field effect device, one uses an effective Bohr radius a* B = h e/(m*e 2 ) 
which includes the dielectric constant e of the medium in which the 
2DEG is created and the effective mass m* of the carriers. Wigner was 
the first to consider the dilute limit where r s becomes large, the Coulomb 
interactions dominate the kinetic energy in determining the wave func- 
tion and the electrons tend to arrange themselves in a regular lattice. 
It may be argued that a lattice configuration is not consistent with the 
translational symmetry characterizing the 2DEG Hamiltonian in the ab- 
sence of a random substrate. This objection can be removed by forming 
a new wave function which is a linear combination of all translations of 
the original lattice: the resultant wave function will have a uniform elec- 
tronic charge density, as symmetry demands, with an unchanged energy. 
Originally, Wigner assumed a bcc electron lattice. It was later shown 
[13] that the hexagonal lattice has a lowest electrostatic energy in two 
dimensions. 

As explained in Refs. [14, 15, 16], when one considers the Hamiltonian 
of N electrons in two dimensions, 

1 " , 2 JL 1 
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r s in 



where the lengths are given in units of l/y/7rn s , the question has been 
from the early days to obtain the asymptotic behaviors of the ground 
state energy Eq around the weak coupling limit (r s « 1): 



£o = ^ + - + 0(lnr s ) (3) 

and around the strong coupling limit (r s >> 1): 

Eo = ^ + 4-2+-2+°( r s 5/2 ), (4) 

r s r -s' r i 

to calculate the coefficients h and /, and to discuss the expected range 
of validity for those asymptotic expansions. Then, one can try to numer- 
ically determine the value of r s where the weak coupling energy exceeds 
the large coupling energy. This can be done at the price of certain ap- 
proximations which are controlled in the two limits and remain more 
uncertain in the middle, and after extrapolating finite size studies to- 
wards the thermodynamic limit. 



The most advanced works in this field are Quantum Monte Carlo 
studies [16, 17] of systems involving many electrons (typically more than 
10 2 ) and assuming two trial wave functions ^t(R) adapted to describe 
either the weak coupling limit, or the strong coupling limit. ^>t(R) ar e 
of the Slater-Jastrow form 

N 

*T(R) = D(R)exp(- J Eu(\r i -r j \)), (5) 

i<j 

where D(R) is a Slater determinant of extended plane waves for weak 
coupling, of localized single-particle orbitals for strong coupling. The 
liquid and crystal pseudopotentials u(r) are repulsive and include in an 
approximate way the effects of electronic correlations. Then a simple 
variational approach, or a more involved fixed node Green's-function 
approach, are used to obtain the energies Eo(r s ) dictated by the cho- 
sen *1>t(R) or by its nodal structure. Both the variational energies and 
the fixed node energies give an upper bound to the exact energy. Com- 
paring for intermediate r s the energies given by the ^t(R) adapted to 
describe the liquid and the crystal, one concludes [17] that there is a 
first order quantum liquid-solid transition at r s ~ 37, with a possible di- 
vision [16, 18] of the liquid phase into a non polarized liquid at small r s 
and a polarized liquid for larger r s . However, the nature of the quantum 
mechanism of melting is still debated, and the possibility of a continuous 
transition has been very recently proposed[19]. The solid is assumed to 
be a frustrated antiferromagnet[20] before becoming ferromagnetic [21] 
at very large r s . The same Monte Carlo method has been used [22] in the 
presence of impurities. The conclusion was that disorder can stabilize 
the solid to weaker values of r s . 

Andreev and Lifshitz have discussed [23] in 1969 the possibility to 
have a more complex intermediate state between the solid and the liq- 
uid, which should be neither a solid nor a liquid. Two kinds of motion 
should be possible in it; one possesses the properties of motion in an 
elastic solid, the second possesses the properties of motion in a liquid. 
This idea comes from a theory of defects in quantum solids. The nature 
of the relevant defects is a complicated issue. Let us give a possible 
example: a vacancy yielded by one electron hopping from the Wigner 
lattice towards some interstitial site. In a classical solid, this defect 
has a certain electrostatic cost and remains localized. In a quantum 
solid, we have in addition the tunneling effect, and if this defect can be 
created, it will be delocalized since the system is invariant under trans- 
lation. Therefore, when r s decreases from the strong coupling limit, the 
increasing band width of the zero point defects of this type may exceed 
their decreasing electrostatic energy cost, leading to two possibilities for 
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intermediate r s : Either the total melting of the solid to directly give a 
liquid, as implicitely assumed for instance in Refs. [19], or a quantum 
floppy solid coexisting with a liquid of delocalized defects, as conjectured 
by Andreev and Lifshitz. A phenomenological FLT theory a la Landau 
of such gapless "delocalized excitations" of a floppy quantum solid has 
been later proposed in Ref. [24]. The discussion of this second possibility 
is one of the central points of this chapter. 

2. DILUTE 2DEG IN FIELD EFFECT 
DEVICES 

A pure 2DEG can be realized by trapping electrons on the surface 
of liquid helium, but it is difficult to reach a sufficient density to study 
the quantum regime. Another possible realization is given by the new 
classes of superconducting cuprates where the electronic motion is es- 
sentially two dimensional. The charge density can be varied by chemical 
doping, and a complex phase diagram is obtained, with insulating, su- 
perconducting and metallic behaviors. 

Eventually, one can create a two dimensional gas of charges (electrons 
or holes) at the interface between two doped semiconductors (GaAs- 
AlGaAs heterostructures), between a semi-conductor and an insulator 
(Si-Mosfet), or very recently [12] between an organic crystal (pentacene, 
tetracene and anthracene) and an insulator. The carrier density can 
be varied by a gate from a very dilute limit towards larger densities. 
In Fig. 1, we mention the remarkable phenomena observed in organic 
(left) or doped semi-conductor (right) field effect devices. A clean inter- 
face may give a high carrier mobility, may allow the observation [26] of 
the fractional quantized Hall effect (FQHE) and may give a measurable 
conductivity in a very dilute limit (typically n s ~ 10 9 — 10 11 carriers 
per cm 2 ) for GaAs heterostructures and organic devices. If the effective 
mass of the carriers is large enough (a condition which is not satisfied by 
the electrons in Ga-As heterostructures) the effective factor r s can be in 
the vicinity of the values where the Fermi- Wigner crossover is expected. 
One of the surprises in those high quality field effect transistors has been 
the observation of a metallic low temperature behavior [4, 12, 27] in a 
certain intermediate range of carrier densities, where a large perpendicu- 
lar magnetic field yields FQHE or magnetically induced Wigner crystals 
[28]. 

3. METAL-INSULATOR TRANSITION 

A recent review of the 2D-MIT can be found in Ref. [4] with an 
extended list of references. We summarize by a few sketches some of the 
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Figure 1. Schematic picture of an organic field effect transistor (upper left) and 
of a GaAs heterostructure (upper right) where the location of the 2DEG or 2DHG 
(symbol plus for the holes) is indicated. The characteristic low temperature behaviors 
are summarized below as a function of the (typical) carrier density n a (cm -2 ): 2dMIT, 
2d metal, FQHE for a sufficient magnetic field, weak localization correction to the 
Boltzmann conductivity, superconductivity (see Ref.[25]). The densities give a typical 
order of magnitude, the observed behaviors depending also on the effective mass of 
the carriers. 



behaviors which have been observed, and which are useful for discussing 
our numerical results. The main surprise was given by the temperature 
dependence of the 2DEG resistivity p(T) around a low critical density 
n c \. As sketched in Fig. 1 (see Fig.l of Ref. [29]), p(T) decreases as 
a function of T when n s < n&, becomes temperature independent at 
n s = n c \ and increases when n s > n c \. A decay is the expected behavior 
for an insulator, while an increase usually characterizes a metal. These 
behaviors occur [30] in a low temperature range 35mk < T < Tp, where 
Tp ~ 0.8— 5-JT are the typical Fermi temperatures of those dilute 2DEGs. 
The temperature increase of p(T) can be large for a 2DEG created in 
a Si-Mosfet (typically one order of magnitude), but remains weak in a 
2DHG created in a GaAs heterostructure. For the densities n s where 
p(T) has a metallic behavior, a parallel field B induces a large positive 
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Figure 2. Resistivity versus temperature without (left) and with (right) a large 
parallel magnetic field B > B sat . n s < n c i (upper curves), n = n c i (middle curves) 
and n > n c i (lower curves). 



magnetoresistance which saturates above a certain field B sa t, as sketched 
in Fig. 3 (see Fig. lofRef. [31] and Fig. 3ofRef. [32]). From small angle 
Shubnikov-de Haas measurements done in a Si-mosfet, it was concluded 
in Ref. [31] that B sa t signals also the onset of full spin polarization. 
Close to the MIT, B sa t is very small and increases as n s — n c \ above 
n c \ [32]. This corresponds to the intermediate values of r s (typically 
3 < r s < 10) where the metallic behavior is observed. When B > B sa t, 
the metallic increase of p(T) disappears, but the I — V characteristics 
sketched in Fig. 4 (Fig 2 of Ref. [29] ) indicates the existence of a critical 
density n C 2 below which a non linearity is observed and above which 
it disappears. The density n c \ and n C 2 are close to each others, if not 
identical when B = 0. The dependence of the characteristic n C 2 as a 
function of a parallel magnetic field B is sketched in Fig. 4. (Fig. 4 of 
Ref. [29]). 

The critical density n c \ does not give a unique critical value for the 
factor r s . Impurity scattering plays a role. For clean systems, one needs 
to have a much larger factor r s than in a dirty system, as sketched in 
Fig. 5. (see inset of Fig. 1 in Ref. [33]). In an (undoped) organic field 
effect transistor, the 2DEG is less scattered by impurities and the MIT 
is seen[27] at an even smaller density (r s ~ 50). 

The surprise caused by this unexpected metallic behavior was mainly 
due to its discrepancy with the scaling theory of localization, which does 
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Figure 3. Resistivity versus parallel magnetic field (left) and saturation field B aat 
as a function of the carrier density (right). 
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Figure 4- Current- voltage non linear characteristics of the insulating phase which 
are suppressed above the MIT (left). Critical density n C 2 above which the non linearity 
of the I-V characteristics disappears as a function of the parallel magnetic field. 



not take into account electron-electron interactions. When the conduc- 
tance g is larger than the conductance quantum e 2 /h, a weak disorder 
perturbative expansion gives for the average conductance a universal log- 
aritmic correction to the Drude conductivity which defavors transport, 
unless there is a sufficient spin orbit scattering. In the weak disorder 
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Figure 5. Critical factor r c s at which the MIT is observed as a function of the inverse 
elastic scattering time. 



limit, one can also take into account the interaction when r s < 1 and 
one obtains additional corrections which reduce transport in a similar 
way. However, the extrapolation of the small r s interaction dependent 
correction to larger r s suggests a possible change of the sign of the cor- 
rections due to the interactions, indicating the possibility of a metallic 
phase in two dimensions, as mentioned by Finkelshtein [34]. Computer 
calculations without interaction and transport measurements at not too 
low densities made in the eighties have confirmed the absence of metallic 
behavior in two dimensions. The difference between the recent experi- 
ments giving a MIT and the former experiments confirming the absence 
of metallic behavior seems to be the quality of the interfaces at which 
the 2DEG is created. This feature makes possible to have a measurable 
conductivity at much lower carrier densities than previously. An inter- 
mediate range of density, where the factor r s is too large (too small) 
to allow expansion in powers of r s (l/r s ) and a weak elastic scattering 
seems to be necessary for observing the metallic behavior. This hypoth- 
esis was supported by Ref. [35] where a study of a 2DHG in a Ga-As 
heterostructure gives a range rf < r s < r^ for having a weak metallic 
behavior in a disordered sample. When r s < rf , one would have weakly 
interacting quasi-particles dominated by Anderson localization when the 
temperature T — ► 0. When r s > rf , one would have a highly correlated 
set of charges. Between rf « 6 and r 



w 



9 in the studied sample, 



10 

a problematic small metallic behavior is observed between two insulat- 
ing behaviors of different nature (Anderson insulator for large densities, 
pinned Wigner solid for low densities). The re-entrant MIT at low r s is 
not easy to observe, if it exists, since the localization length of a clean 
device can be very large, and the observation of a possible re-entrant 
insulating behavior at high densities can require very low temperatures. 
This is why many observations of a MIT have been reported for r s ~ 10, 
while very few experiments give a possible re-entrant MIT at r s ~ 3. 
Moreover, more recent works [5, 7, 9]) put doubts about the reality of 
this intermediate metallic behavior when the temperature goes to zero, 
since the effect of a weak perpendicular magnetic field can be described 
by usual weak localization theories, even for values of r s as large as 15 
(see Fig. 6 taken from Ref. [5]). The hypothesis of a certain temper- 
ature dependent screening was suggested for explaining the anomalous 
temperature dependence, and it was proposed that usual quantum in- 
terferences should drive at possibly very low temperatures the system to 
the formerly expected insulating behavior. Recent measures performed 
down to 5mK do not confirm [36] this hypothesis. The estimate of the 
phase breaking length L^ which is traditionally done for estimating the 
low temperature dependence of the resistance from a zero temperature 
theory leads again to the famous problem of the saturation [37] of L^ 
when T — ► 0, problem leading also to many possible and controversial 
explanations. 
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Figure 6. Usual weak localization behavior of the conductance induced by a per- 
pendicular magnetic field in a 2DHG created in Ga-As heterostructure. 
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In summary, a significant metallic behavior can be seen using a 2DEG 
created in a Si-mosfet, while a weaker one occurs in a 2DHG created 
in a Si-Ge quantum well or a Ga-As hetrostructure. Nevertheless, in 
the latter system, the study of the compressibility gives complementary 
signatures [38, 39, 40] of a possible quantum phase transition. Local 
compressibility measurements show that the system is more homogenous 
in the intermediate metallic phase than in the low density insulating 
phase. Very recently, the possibility that the MIT would be accompanied 
by a magnetic transition has been suggested [32, 41]. 

4. LATTICE MODEL 

The previous experimental observations lead us to numerically revisit 
the Fermi- Wigner crossover using a two dimensional model describing 
N particles on L x L square lattice with periodic boundary conditions 
(BCs), i.e. with a torus topology. The most general Hamiltonian TC of 
the lattice model we will focus on reads, 

i,cr i' 

^E^ + ^E™.,^ (6) 

iyti' 

where the operators Cj iCr (c ia ) destroy (create) an electron of spin a at 
the site i and nj iCT = c ia Ci :ff . 7i consists of 

■ a hopping term —t that couples nearest-neighbor sites, and ac- 
counts for the quantum kinetic energy, 

■ the pairwise electron-electron interaction, which itself consists of 
a 2U Hubbard repulsion when two electrons are at the same site 
i with opposite spins and a U/\i — i'\ spin independent Coulomb 
repulsion when they are separated by a distance \i — i'\ (smallest 
distance between the sites i and i' on a square lattice with periodic 
BCs), 

■ on site random potentials V{ which are uniformly distributed inside 
the interval [-W/2,W/2]. 

The clean system is obtained when the disorder strength W is set to zero. 
In our model with a lattice spacing a, h /(2m*a 2 ) — ► t, e 2 /(ea) — ► U, 
such that the factor r s becomes: 

(7) 



y/irn s a* B 2t^/im e ' 
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Figure 7. Structure of the Hamiltonian matrix for 4 electrons in a 6x6 square lattice. 
The size of the different diagonal non zero sub-blocks are N(S Z ) — 396900, 257040 
and 58905 for S z = 0, ±1, ±2 respectively. 



for a filling factor n e = N/L 2 . This dimensionless ratio r s will allow 
us to compare our results obtained as a function of U for a fixed filling 
factor n e and the experimental results obtained as a function of n s . 

We denote S and S z the total spin and its component along an ar- 
bitrary direction z. Since [S 2 ,7i] = [S Zl 7i] = 0, TL can be written in 
a block-diagonal form, with N + 1 blocks where S z = —N/2, . . . , N/2 
respectively. When B = 0, there is no preferential direction and the 
groundstate energy Eq does not depend on S z . For a groundstate of to- 
tal spin S, TL has 25+1 blocks with the same lowest eigenenergy Eq(S 2 ) 



since E (S 2 ) = E (S 2 ,S Z 



-S,-S+i,...,S-i,S. Therefore, the number 



Nt> of blocks of different S z and of same lowest energy gives the total 
spin S = (Ni> — l)/2 of the groundstate. 

If N and L are small enough, the ground state and the first excitations 
can be exactly calculated using Lanczos algorithm. Otherwise, certain 
approximations are unavoidable. Let us focus on the case N = 4 and 
L = 6 where the structure of TC is given in Fig. 7 for different S z . 
Without magnetic field, we have the symmetry ±S Z , and we have only 
to diagonalize the three sub-blocks with S z > 0. TC(S Z = 2) corresponds 
to fully polarized electrons (spinless fermions) where the orbital part of 
the wave-functions is totally anti-symmetric. 
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From exact diagonalization for small systems and using approxima- 
tions for larger systems and weak coupling, we will study: 

■ the lowest eigenenergies E n (S z ) of the S z sub-blocks and the cor- 
responding eigenvectors \^ n (S z ) >■ n = 0, 1, 2, . . . corresponds to 
the states ordered by increasing energies. 

■ The local persistent currents J(i) created at a site i by an Ahronov 
Bohm flux <j) which is enclosed along the longitunal /-direction 
as sketched in Fig. 8. The flux </> can be included by taking 




Figure 8. 2D Torus with N electrons enclosing an Aharonov-Bohm flux <f> — BL% 

appropriate longitudinal BCs (antiperiodic BCs corresponding to 
= 7r in our convention). The BCs along the transverse ^-direction 
remain periodic. The J{i) are vectors defined by their longitudinal 
and transverse components (J(i) = (Ji,i, Ji,t)), or by their angles 
0i = arctan(Jj )i /Jj^) and their absolute values Jj = \Ji\. The 
longitudinal component Jn of an eigenstate |^) is defined as 



Ji,i = 2Im(*|cJ r+liiy Ci Bi i s exp(i0/L)|tt) 
and Ji t t is given by 

J i , t = 2Im<¥|4 jVh iCi.,ij¥>. 



(8) 
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The total current P n > of the n many-body wavefunction |\l/ n ) of 
energy E n has a total longitudinal component Ii(n) given by 



iiinm = -f? 



J2iJi,l(n) 



=~4> L 

which will be calculated for cp = ir/2. 

The crystallization parameter 7 defined using the function C(r) = 
N^ 1 J2i PiPi-n where pi = (^\ni\^) is the electronic density of the 
state |\&) at the site i. The crystallization parameter 7 is given by 

7 = max C(r) — min C(r) (11) 

Note that 7 = 1 when the N particles are localized on N lattice 
sites and form a rigid solid and when they are extended on the 
L 2 sites and form an homogenous liquid. 

The participation ratio % = N 2 (J^ pf)^ 1 , which gives the typical 
number of lattice sites occupied by an eigenstate |\& >. 

The spectral parameter r\ which characterizes the level repulsion. 
Uncorrelated spectra exhibit Poisson statistics. Correlated spec- 
tra can be described by Random Matrix Theory with Wigner- 
Dyson (W-D) statistics. For the one body spectra, the distribution 
P(s) of the normalized energy spacings between consecutive levels 
has two different forms when L — > 00: the Poisson distribution 
Pp(s) = exp(— s) if the wavefunctions are localized, the Wigner 
surmises P{^(s) = (tts/2) exp(— its 2 /A) with time reversal symme- 
try (TRS) and P^{s) = (32s 2 /vr 2 ) exp(-4s 2 /vr) without TRS, if 
the wave functions are extended. For the normalized iV-body en- 
ergy spacings s n = (E n+ i — E n )/ < E n+ \ — E n > (the brackets 
denote ensemble average), we define a spectral parameter: 



var(P( g ))-var(P^ ) ( g )) 
v a T(Pp(s))-v a r(P^ U \s)) 



^o,u) = — — — ,%m, . ' ( 12 ) 



(rjo with TRS, r/jj in the absence of TRS, for instance when <f> = 
7r/2)). var(P(s)) denotes the variance of P(s). The spectral pa- 
rameter rj = 1 when P(s) = Pp(s) and 7/(0 [/) = when P(s) = 

P& U \s). 

The Zeeman energy necessary to polarize a non magnetized cluster. 
A parallel magnetic field B does not induce orbital or Aharonov- 
Bohm effects, but defines the z-direction and removes the S z de- 
generacy by the Zeeman energy —gpBS z . The ground state energy 
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and its magnetization are given by the minimum of Eq(S 2 , S z , B = 
0) — g(iBS z . For a S = groundstate without field, the value B* 
for which Eq(S z ) — gfiB*S z = Eq(S z = 0) defines the field neces- 
sary to polarize the system to S > S z . If one studies N = 4 
electrons, the total Q2 = Eq(S z = 2) — Eq(S z = 0) and partial 
Qi = Eq(S z = 1) — Eq(S z = 0) polarization energies give the Zee- 
man energies necessary to yield S = 2 and S = 1 respectively for 
a cluster with S = 0. 

6. INTERMEDIATE COUPLING REGIME 

FOR SPINLESS FERMIONS AND WEAK 
DISORDER 

We first consider an ensemble of disordered clusters with L = 6 and 
N = 4. The ground state (GS) and the first excitations of the fully po- 
larized sub-block (S z = 2, spinless fermions) of the Hamiltonian matrix 
shown in Fig. 7 have been obtained using the Lanczos algorithm. The 
statistical ensemble typically includes 10 3 — 10 4 samples obtained from a 
disorder distribution with W = 5. This is a relatively weak disorder for 
which one has quantum diffusion when r s = inside the small clusters 
(no Anderson localization). 

6.1. GROUND STATE 

We summarize in this subsection the main results published in Ref. 
[42] and complementary unpublished results. As one switches on U, 
a first characteristic threshold rf can be identified by looking at the 
average total longitudinal persistent current Ii of the GS at <f> = tt/2, 
and comparing the exact quantity with the Hartree-Fock (HF) approx- 
imation (see appendix). Below rf ~ 5, the mean field approximation 
reproduces the exact Ii, but strongly underestimates Ii above rf. This 
sharp breakdown of the HF approximation shown in Fig. 9 means that 
strong correlation effects occur above rf , such that the shift the GS en- 
ergy when the BCs are changed cannot be obtained assuming the best 
possible SD for the ground state. 

A closer investigation of the persistent currents on a typical sample 
gives three regimes, as shown in Fig. 10. See also Ref. [43]. For weak 
coupling, the local currents flow randomly inside the cluster, due to 
elastic scattering on the site potentials. For intermediate coupling, the 
pattern of the persistent currents becomes oriented along the shortest 
direction enclosing d). For large coupling, the oriented currents vanish. 
Ref. [44] gives a detailed study of the large coupling limit where one can 
use perturbation theory for having the sign and the magnitude of Ii. 
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Figure 9. Ensemble average longitudinal GS current < Ii > as a function of r s 
for N = 4, L — 6 and W = 5. Exact values (filled symbols) and HF values (empty 
symbols) . 
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Figure 10. Map of the local persistent currents in a given sample for small (left), 
intermediate (center) and large (right) values of r a . 



If one looks at the distribution of the angles 0i of the local currents, 
one can see in Fig. 11 that the currents are randomly scattered without 
interaction, and that they become aligned when one goes to the strong 
coupling limit. The ensemble average value < \9\ > allows us to quantify 
the progressive change. If p(6) = l/(2ir), < \8\ >= n/2, a value obtained 
for the low ratios r s . At large r s , < \9\ >^ 0. The ratio r s at which 
the local currents cease to be oriented at random is consistent with the 
critical ratio rf where the HF approximation breaks down. 
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Figure 11. Left: Distribution Pe of the GS local current angles for r s — (dotted 
line), r s — 6.3 (full line) and r s = 42 (dashed line). Right: Ensemble average angle 
< \8\ > as a function of r s . 



By studying the average amplitude of the local currents, one can see 
in Fig. 12 that < J« > is essentially independent of r s up to a second 
threshold r^ ~ 10 which exceeds rf. Moreover, comparing in Fig. 12 
the GS average crystallization parameter < 7 > and < Jj >, on can 
see that the suppression of the persistent currents coincides with the 
formation of a solid Wigner molecule inside the disordered clusters. 
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Figure 12. Averages of the GS crystallization parameter < 1 — 7 > (left scale, 
empty symbols) and GS local current amplitude < Ji > (right scale, filled symbols) 
as a function of r s . 
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The response of the ground state to an enclosed Aharonov-Bohm flux 
shows us that an intermediate correlated regime takes place between the 
Fermi limit and the Wigner limit, when typically 5 < r s < 10. 



6.2. 



LOW ENERGY EXCITATIONS 



In Ref. [45], the low energy excitations of the same clusters have been 
studied, notably their statistics when the microscopic configurations of 
the random substrate are changed. For intermediate ratios r s , the GS 
and the 8 first following low energy excitations are characterized by ori- 
ented non random persistent currents and do not exhibit Wigner-Dyson 
(W-D) spectral statistics. Above those states, when the excitation en- 
ergy e exceeds an energy of the order of the Fermi energy ef, the local 
currents become randomly oriented and the levels obey W-D statistics. 
Incidentally, let us note that the metallic behavior observed for interme- 
diate couplings disappears also when the temperature exceeds a temper- 
ature of the order of the Fermi temperature. 
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Figure 13. Left: Spacing distribution P(s) for r s = 6.3 when cj> — (o) and <j> = 7r/2 
(•) for energy levels above the 9 first levels (excitation energies 1.4 < e/ef < 1-9), 
compared to Pq (s) (dashed line) and P]f (s) (continuous line). Right: Spectral 
parameter r\u as a function of r a for the first excitation so (circles), S2 — S4 (squares) 
and sio — S20 (diamonds). 



In Fig. 13 (left) one can see that the spacing distribution P(s) calcu- 
lated for r s = 6.3 using the low energy levels except the 9 first levels (n = 
10, ... , 20), is given by the Wigner surmise, with an orthogonal-unitary 
crossover when one turns on an Aharonov-Bohm flux <fi = — ► ir/2. This 
corresponds to excitation energies 1.4 < e/ef < 1-9. Taking <\> = it/2, 
the variation of the spectral parameter rju as a function of r s is given for 
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the successive level spacings in Fig. 13. The first excitation is described 
by the Wigner surmise (rjjj = 0) without interaction but becomes more 
and more Poissonian when r s increases. The spacings characterizing the 
levels with n = 10, . . . , 20 have an opposite behavior. For r s = 0, those 
excitations being the sum of more than one single-electron excitation 
are essentially uncorrelated, but become correlated for intermediate r s 
(rju ~ 0) before being again uncorrelated at larger r s . 

In summary, when one considers the low energy spectral statistics, a 
complementary signature of an intermediate regime is obtained, given by 
W-D statistics and randomly oriented local persistent currents outside 
the 9 first states for which the absence of W-D statistics for intermediate 
r s is accompanied by a non random orientation of the persistent current 
angles 9 (see Ref. [45]). This behavior does not appear for weak and 
strong couplings, where the low energy spectral correlations decrease as 
the excitation energy increases. 

6.3. INTERMEDIATE LIQUID-SOLID 

REGIME IN THE CLEAN LIMIT 

The fact that the 9 first states do not display quantum ergodicity 
for intermediate coupling and weak disorder suggests the existence of 9 
low energy collective excitations. A collective motion cannot be due to 
impurity scattering and should come from the corresponding clean limit. 
This limit has been investigated in Ref. [46]. 

When W = 0, one has a single system which remains invariant under 
rotation of angle tt/2 and under translations and reflections along the 
longitudinal x and transverse y directions. Invariance under translations 
implies that the momentum K is a good quantum number which remains 
unchanged when U varies. The symmetries imply that the states are 
fourfold degenerate if K ^ and can be non degenerate if K = 0. 

When C/ = 0, the states are Nh plane wave Slater determinants (SDs) 

^1(4)^1(3)^1(2)^1(1) I ^ ■ > ' wnere ^k(n) crea t es a particle in a state of mo- 
mentum k(p) = 2ir(p x ,Py)/L (p x ,y = I, • • • j L) and |0 > is the vacuum 
state. For N = 4 and L = 6, Njj = 58905. The low energy eigenstates 
are given by the following plane wave SDs: 

■ 4 degenerate ground states (GSs) \Kq((3) > {(5 = 1, . . . , 4) of energy 
E Q (U = 0) = —13* and of momenta K = (0, ±7r/3) and (±vr/3,0). 

■ 25 first excitations of energy E\(U = 0) = —12* out of which 4 
plane wave SDs \K\{(3) > will play a particular role for describing 
the intermediate GS. They correspond to a particle at an energy 
—At with k{l) = (0, 0), two particles at an energy —3* and a fourth 
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particle of energy — It with momenta such that J2j=2 Mi) = 0- One 
has k(2) = (0,±7r/3), fc(3) = (±tt/3,0) and fc(4) = (=Fvr/3, =Fvr/3) 
or fc(2) = (0,=ftt/3), fc(3) = (±vr/3,0) and fc(4) = (=Ftt/3, ±tt/3). 

■ 64 second excitations ^2(0;) > of energy E^{U = 0) = —lit. 

■ 180 third excitations |K3(a) > of energy ^(C/ = 0) = —104. 

■ 384 fourth excitations of energy E^(U = 0) = —94 out of which 16 
plane waves SDs \K^(5) > will play a particular role for describing 
the intermediate GS. They are given by the condition that the to- 
tal momentum is zero, which selects 20 SDs out of which 4 where 
the single particle state of energy —34 is not occupied do not con- 
tribute. The \K±(8) > are 16 SDs of energy —94, given by 8 SDs 
where the particles have energies — 44, — 34, —24, 04 respectively and 
by 8 other SDs where the particles have energies — 34, —34, — 2t, —t 
respectively. 

When t = 0, the states are Nh Slater determinants qric[cj|0 > built 
out from the site orbitals. The configurations ijkl correspond to the 
Nh different patterns characterizing 4 different sites of the 6x6 square 
lattice. The low energy part of the spectrum is made of the following 
site SDs: 

■ 9 squares |Sb(7) > (I = 1, - . . , 9) of side a = 3 and of energy 
E (t = 0) w 1.80*7. 



36 parallelograms \Si(I) > of sides (3, v 10) and of energy ~ 1.8577. 



■ 36 other parallelograms |5*2(7) > of sides (vlO, v 10) and of energy 
w 1.9717. 

■ 144 deformed squares 15*3(7) > obtained by moving a single site of 
a square \Sj > by one lattice spacing and of energy w 2U. 

For the first low energy states, the crossover from the U = eigenbasis 
towards the t = eigenbasis is shown in Fig. 14 when one increases the 
ratio r s . If we follow the 4 GSs Eo(r s = 0) (Kq 7^ 0), one can see a 
first level crossing at rf sa 9.3 with a non degenerate state (Kq = 0) 
which becomes the GS above rf, followed by two other crossings with 
two other sets of 4 states with Ki / 0. When r s is large, 9 states 
coming from E\ (r s = 0) have a smaller energy than the 4 states coming 
from Eo(r s = 0). The degeneracies ordered by increasing energy become 
(1, 4, 4, 4, . . .) instead of (4, 25, 64, . . .) for r s = 0. Since the degeneracies 
are (9, 36, 36, . . .) when t = 0, these 9 states give the 9 square molecules 
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Figure 14- As a function of r s , low energy part of the spectrum exhibiting a GS 
level crossing at rf . Inset: two first level spacings Ai/i (dashed) and A2/t (dotted) 
which become equal at r^ and the perturbative result Ai/i = A2/t 
when r s — > oo (dot-dashed). 
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|5o(/) > when r s — ► oo. When r^T 1 is very small, the first 9 states 
correspond to a single massive molecule free to move on a restricted 
3x3 lattice, the single non frozen degree of freedom in this limit being 
the location Rj of the center of mass of the |So(-f) >• One has an 
effective hopping term T oc £r,T 3 when N = 4 and the total momentum 
is quantized (K[(I) = 2irpi/3 and K t {I) = 27rpt/3 being its longitudinal 
and transverse components respectively with pit = 1, 2, 3). For a square 
lattice at a filling factor 1/9, the Rj are indeed located on a periodic 3x3 
square lattice. This is an important simplification of our model. This 
gives 9 states of kinetic energies given by — 2T (cos K~i(I) + cos K t (I)). 
The kinetic part of the low energy spectrum is then — 4T, — T, +2T with 
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degeneracies 1,4,4 respectively. This structure with two equal energy 
spacings Ai and A2 appears (inset of Fig. 14) when r s is larger than the 
crystallization threshold rf ~ 28. Above r^ , to create a defect in the 
rigid molecule costs a high energy available in the 10 excitation only. 
We have seen in the previous section that the 9 first levels do not obey 
Wigner-Dyson statistics at intermediate r s when a random potential is 
added, in contrast to the following levels. The study of the clean limit 
gives us the explanation. The two characteristic thresholds rf (level 
crossing) and r^ (9 first states having the structure of the spectrum of 
a single massive molecule free to move on a 3 x 3 square lattice) can also 
be detected by other methods given in Ref. [46]. 




Figure 15. Ground state projections Po(r 3 ) onto a few plane wave SDs, given by 
the A\K (P) > (empty circle), the 4 \Ki(/3) > (empty square), the 64 \Ki(a) > (filled 
diamond), the 180 \K2(a) > (x), the 16 \Ki{8) > (asterisk) respectively, as a function 
of r 3 . 



To understand further the nature of the intermediate GS, we have 
projected the GS wave functions | v &o( r s) > over the low energy eigen- 
vectors of the two eigenbases valid for U/t = (Fig. 15) and for t/U = 
(Fig. 17) respectively. 
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Let us begin by studying the GS structure in the U = eigenbasis. 
Below rf , each of the 4 GSs |\l/g(r s ) > with Kq ^ has still a large 
projection 

4 

P (r s ,0) = ]T | < V%(r a )\K (j3) > | 2 (13) 

p=i 

over the 4 non interacting GSs. There is no projection over the 25 first 
excitations and smaller projections Pq{t s ,2) and Pg( r s>3) over the 64 
second and 180 third excitations of the non interacting system. Above 
rf, the non degenerate GS with Kq = has a large projection 

4 

P (r s , 1) = E I < *o(r.)|tfi(/3) > | 2 (14) 

/3=1 

which is equally distributed over the 4 excitations | K\ (/3) > of momen- 
tum il"i = and a second significant contribution 

16 

P (r s ,4) = Y / \<^o(r.s)\K 4 (5)>\ 2 (15) 

5=1 

given by its projection onto the 16 plane wave SDs \K^(5) > belonging 
to the fourth excitation of the non interacting system. Above rf, its 
projections onto the 4 \K$(0) >, the 21 other first excitations and the 
second and third excitations of the non interacting system are zero or 
extremely negligible. 

The total GS projection 

P*(r s )=P (r s ,l)+P (r s ,4) (16) 



onto the 4 \Ki(/3) > and 16 \K4(5) > is given in Fig. 16 when r s > rf . 
This shows us that a large part of the system remains an excited liquid 
above rf, given by a special rule of occupation of the one particle plane 
wave states. The occupation of the one body states is very different from 
the usual Pauli rule after the level crossing rf. A necessary, but non 
sufficient condition for a plane wave SD to significantly contribute to the 
zero momentum GS is of course to have a zero total momentum. Those 
projections decrease as r s increases and become negligible in the strong 
coupling limit. A complete GS description in this limit will require more 
and more plane wave SDs. 

We now study the GS projections P^ onto the t = eigenbasis. The 
GS projection 

9 

Poo(r s , 0) = J2\< n(r s )\So(I) > | 2 (17) 
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Figure 16. Ground state projection P (r s ) (asterisk) and Ptx>(r s ) (empty square) 
onto the subspace spanned by the low energy plane wave and site SDs respectively, 
and total GS projection P (filled circle) onto the re-orthonormalized basis using the 
low energy eigenvectors of the two limiting bases. 



onto the 9 square site SDs |<5o(-0 > is given in Fig. 17, together with 
the GS projection Poo{r s , J) onto the site SDs corresponding to the J th 
degenerate low energy excitations of the t = system. The total GS 
projection 



P l oo(r s 



p=0 



,p) 



(18) 



onto the 9 squares | So (I) > , the 36 parallelograms | S± (I) > , the 36 other 
parallelograms IS2CO > an d the 144 deformed squares (^(J) > is given 
in Fig. 16 when r s > rf. This shows us that the ground state begins 
to be a floppy solid also above rf . When r s increases, Poc{r s , 0) goes to 
one, and the ground state is a simple rigid square Wigner molecule. 

The site SDs and plane wave SDs are not orthonormal. After re- 
orthonormalization, the total projection P of |^o( r s) > over the sub- 
space spanned by the 4 \K\{(5) > and 16 |i^4(<5) > and 225 site SDs 
of lower electrostatic energies (9 squares, 36 + 36 parallelograms, 144 
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Figure 17. Ground state projection Poo(r s ) onto a few site SDs, given by the 9 
squares |So(/) > (filled square), the 36 parallelograms Si (J) > (asterisk), the 36 
parallelograms |S I 2(/) > (diamond), and the 144 deformed squares \Sz{I) > (left 
triangle) respectively, as a function of r 3 . 



deformed squares) are given in Fig. 16. One can see that |^o( r «) > 
is almost entirely located inside this very small part of a huge Hilbert 
space for intermediate r s , spanned by no more than 245 low energy SDs 
of different nature, and adapted to describe a solid entangled with an 
excited liquid. 

From the study of the GS projections emerges the conclusion that a 
minimal description of the intermediate GS requires to combine the low 
energy states of the two limiting eigenbases. In this sense, the GS is 
neither solid, nor liquid, but rather the quantum superposition of those 
two states of matter. This is strongly reminiscent of the conjecture 
proposed in Ref. [23] for the quantum melting of solid Helium in three 
dimensions. It suggests possible improvements of the trial GS to use 
for intermediate r s in variational or path integral quantum Monte Carlo 
approaches. Instead of using Jastrow wave functions improving the plane 
wave SDs for the liquid or the site SDs for the solid, or their nodal 
surfaces, it will be interesting to study if a combination of the two, with a 



26 

more complex nodal structure, and describing a solid-liquid regime does 
not lower the GS energy for intermediate r s . A positive answer would 
confirm that an unnoticed intermediate solid-liquid phase does exist in 
the thermodynamic limit for fermionic systems in two dimensions. 

6.4. MAGNETIC SIGNATURE OF THE 

INTERMEDIATE REGIME FOR WEAK 
DISORDER 

We return to the study of weakly disordered samples when the spin 
degrees of freedom are included. Their role and the consequences of 
an applied parallel magnetic field which align only the spins without 
inducing orbital effects, have been the subject of Ref.[47]. A statistical 
ensemble of matrices having the structure given in Fig. 7 have been stud- 
ied for W = 5, N = 4 and L = 6 providing complementary signatures 
of a particular intermediate behavior. 

When r s = 0, the two one body states of lowest energy are doubly 
occupied and S = (S = 1/2 if N is odd). To polarize the 5 = ground 
state to S = 1 corresponds to the transition of one electron at the Fermi 
energy and costs an energy equal to the one body level spacing. p(Qi) 
is then given by the spacing distribution p(s) between consecutive one 
body levels, the Wigner surmise P£(s) in the diffusive regime. When r s 
is large, the 4 electrons occupy the four sites Cj j = 1, . . . , 4 of the square 
configuration |<So(-0 > °f side a = 3 with the lowest substrate energy 
J2j=i v cj- The ground state in this limit becomes \^ c >= Ilj=i c t ,<r |0 > 
with a spin independent energy E c . This square can support 2^ = 16 
spin configurations. We summarize the main results of a perturbative 
expansion of E c in powers of t/U. The spin degeneracy of E c is re- 
moved by terms of order £(£/[/) 2a_1 , which is the smallest order where 
the 16 spin configurations can be coupled via intermediate configura- 
tions allowing a double occupancy of the same site. Therefore, 2a — 1 
is the order where the perturbation begins to depend on S z and Q\ as 
Q2 oc t(t/U) 2a ~ l — ► when t/U — ► (we have numerically checked this 
decay when r s > 100). Moreover, the correction to E c depending on S z 
and oc t(t/U) 2a ~ l is given by an effective antiferromagnetic Heisenberg 
Hamiltonian. The S = ground state for large r s correspond to 4 elec- 
trons forming an antiferromagnetic square Wigner molecule. However 
Qi and Q2 are very small when r s is large, and the antiferromagnetic be- 
havior can be an artefact due to the square lattice. Without impurities 
and in a continuous limit, a quasi-classical WKB expansion [21] valid for 
very large values of r s shows that 3 particle exchanges dominate, leading 
to ferromagnetism. Recent Monte-Carlo calculations [20] suggest that 
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Figure 18. Distributions p of the polarization energies Q\ and Q2 at different values 
of r s . Upper left: p(qi) at r 3 — (circle) and 2.5 (square) where q\ = Qi/ < Q\ >. 
The continuous line is the Wigner surmise. p(— log Qi) (filled circle) and p(— log Q2) 
(empty diamond) at r s — 2.5 (upper right, right scale) 5.8 (lower left, left scale) and 
16.8 (lower right, right scale) respectively. The thick bars (put at right edge of the 
figures) give the peaks 5(Qi) of the bimodal p(Qi). The continuous lines are normal 
fits. 
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the crystal in the continuous limit becomes a frustrated antiferromagnet 
closer to the melting point. 

The perturbative corrections oc t(t/U) 2a ~ 1 depend on the random 
variables Vi via Ilj^i {E c — Ej)~ l where the Ej are the classical energies 
of the intermediate configurations. Ej is the sum of an electrostatic 
energy and of a random substrate energy E S (J) = J2k=i v J(k)- Due to 
the high order 2a— 1 of the correction, a normal distribution for Ej leads 
to a log-normal distribution for IIj=i (E c — Ej)~ l . Therefore p(Ai) and 
p(A2) should be log- normal when r s is large. 

p(Qi) is given in Fig. 18 for different r s . The expected Wigner sur- 
mise takes place for r s = 0. A small interaction quickly drives p(Qi) 
towards a bimodal distribution, with a delta peak at Q± = and a main 
peak centered around a non zero value of Q\. The delta peak gives the 
probability to have spontaneously magnetized clusters with 5 = 1. The 
main peak gives the field B necessary to create 5 = 1 in a cluster with 
5 = 0. The logarithmic scale used in Fig. 18 (upper right) underlines 
the bimodal character of the distribution and confirms that the main 
peak becomes log-normal when r s is large. The distribution of Q2 is not 
bimodal: a fully polarized cluster has never been seen when B = 0. Qi 
becomes also log-normally distributed when r s is large. 

In Fig. 19, the fraction M of clusters with 5 = 1 at B = is given 
as a function of r s . One can see the mesoscopic Stoner instability (see 
appendix) taking place at r s ~ 0.35. The Stoner mechanism should 
eventually give fully polarized electrons. This is not the case, the in- 
crease of M breaks down when r s = rf ~ 2.2, a value where the 
Stoner mechanism and hence the HF approximation break down. In the 
same clusters, we have seen that the HF approximation fails to describe 
the persistent currents of the fully polarized sub-block of TC (spinless 
fermions) when r s > rf ~ 5. rf takes a smaller value rf when the 
spin degrees of freedom are included. Above rf s , M regularly decreases 
to reach a zero value for r^ s ~ 9 where an antiferromagnetic square 
molecule is formed. In the intermediate regime, there is a competition 
between the Stoner ferromagnetism and the Wigner antiferromagnetism. 
Since the 5 = clusters are characterized by log-normal distributions, 
the ensemble averages < log Q\ > and < log Q2 > (without taking into 
account the 5 = 1 spontaneously magnetized clusters) define the typical 
fields B necessary to yield 5 = lor5 = 2ina5 = cluster. Fig. 
19 provides two magnetic signatures confirming the existence of a novel 
intermediate regime between the Fermi glass (r s < rf s ) and the Wigner 
glass (r s > rf s )'- < logQi > becomes roughly independent of r s , while 
Q 2 oc r~ 2 . 
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Figure 19. As a function of r 3 , fraction M of clusters with S = 1 at B — (filled 
diamond, right scale), partial < logQi > (filled circle, left scale) and total < logQ2 > 
(empty triangle, left scale) energies required to polarize S — clusters to S — 1 and 
S — 2 respectively. The straight line corresponds to 0.25 — 21ogr s . 
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Figure 20. As a function of r s , ensemble averages of the crystallization parameter 
< l{S z ) > (right) and of the numbers of occupied sites < x{Sz) > (left). S z — 
(circle), S z — 1 (square) and S z — 2 (triangle). The arrows indicate the effect of a 
parallel magnetic field. 



The number of occupied sites x(Sz) depends on S z for small r s and be- 
comes independent of S z for large r s . At r s = 0, ^q{S z = 2, 1, 0) occupy 
respectively 4, 3, 2 one body states while the Wigner molecule occupies 
4 sites only at large r s . The ensemble average < x{Sz = 2) > shown in 
Fig. 20 is maximum when r s = and decays as r s increases, suggesting 
the absence of derealization for the polarized system. The non polarized 
system behaves differently, since < x(-S , 2 = 1, 0) > first increase to reach 
a maximum ~< x{$z = 2,r s = 0) > before decreasing. In Fig. 20 one 
can see also from the curves j(r s ) that charge crystallization is easier 
when the clusters are polarized than otherwise. The arrows indicated in 
Fig. 20 underline two consequences of a parallel field B: smaller number 
of occupied sites and smaller crystallization threshold. 



6.5. 



NUMERICS VERSUS EXPERIMENTS 



One can question whether the comparison between numerical simula- 
tions based on a small lattice model with a few particles and measures 
using 10 11 electrons per cm? makes sense. It is nevertheless interesting 
to point out certains analogies in the obtained thresholds and behaviors. 



The ratios rf where the rigid Wigner molecule is formed in small 
clusters are of the order of the ratios rj where one observes the 
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MIT. One finds a value ~ 10 for a weakly disordered cluster which 
increases to larger values when one goes to the clean limit (ry ~ 
30 — 40). The r^ numerically obtained for larger disorder (W = 
5, 10, 15) in Ref.[42] are roughly independent of W when W > 5 
and reproduce the flat part of the experimental curve sketched in 
Fig. 5. 

■ The large magnetoresistance sketched in Fig. 3 is consistent with 
the reduction of \ yielded by a parallel magnetic field as indicated 
in Fig. 20 (right). 

■ The shift of the critical threshold rf when one polarizes the elec- 
trons with a parallel magnetic field (Fig. 20 - left) is consistent 
with the shift of the critical densities sketched in Fig. 4. 

■ A polarization energy Q2 oc r~ 2 is consistent with the n s depen- 
dence of B sa t given in Fig. 3 for similar values 3 < r s < 9 [29] . 

■ A threshold at rf ~ 3 — 5 given by the study of the persistent 
currents and of the magnetization corresponds to the density at 
which the metallic behavior ceases to occur, according to Ref. [35]. 

■ It can be argued that the usual small negative magnetoresistance 
yielded by a perpendicular magnetic field for intermediate r s (Fig. 
6) would mean that transport is due to the gapless excitations of 
the intermediate floppy quantum solid. The propagation of such 
excitations could be reduced by the usual weak localization cor- 
rections when one has elastic scattering by impurities. 

6.6. MESOSCOPIC 2D ELECTRON SYSTEMS 

CONFINED IN HARMONIC TRAPS 

For the important issue of the melting of a macroscopic 2D Wigner 
glass (crystal without disorder) and a possibly associated 2D-MIT, the 
study of finite size systems can give hints, but requires to be comple- 
mented by a systematic study of the scale dependence of the finite size 
effects at a given electron density. Before reviewing some results where 
such a finite size scaling analysis has been done, let us underline that the 
formation of the mesoscopic Wigner molecule in a few electron system is 
by itself an important issue. One can use a few electrons [48] confined in 
a quantum dot or a few ions [49] trapped by electric and magnetic fields. 
Increasing the size of the trap, a crossover from independent-particle 
towards collective motion can be observed. Moreover, a localization- 
delocalization transition has been observed [50] in a quantum dot using 
single-electron capacitance spectroscopy and increasing the number of 
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trapped electrons. Considering electrons confined in a harmonic trap, 
and not on a 2D torus, other numerical studies have reached the same 
crucial conclusion, namely that between the Fermi system and the "clas- 
sical" rigid molecule, there is an intermediate regime of a floppy weakly 
formed Wigner molecule. This was shown for instance in an exact study 
of a two-electron artificial atom [51] and in a Monte Carlo study [52] 
of a few electron system. In Ref. [52], this new regime corresponds to 
particles having a radial ordering on shells without correlated intershell 
rotation. This was attributed to the special symmetry of the harmonic 
trap and to the imposed density gradient. Our studies show that meso- 
scopic Wigner crystallization takes also place in two stages when the 
particles are confined on a 2D torus, with a uniform electron density. 

7. QUANTUM PHASE TRANSITION FOR 

WEAK COUPLING AND STRONG 
DISORDER 

In this Section we consider the case of spinless fermions for strong dis- 
order and low filling factors N/L 2 . In contrast to the previously studied 
systems, the relevant length scale without interaction, the one-body lo- 
calization length Li, is not only smaller than the system size, but also 

— 1/2 

smaller than the average distance between particles, a oc n s . Trun- 
cated bases built out from single-particle wavefunctions or Hartree-Fock 
orbitals will be used to study larger systems. These approximations re- 
main reliable when the interaction strenght U does not exceed the hop- 
ping term t, but become uncontrolled in the strongly correlated regime 
U 3> t. They allow us to study the finite size effects in the limit of rela- 
tively weak coupling. They show the existence of a second order phase 
transition between the weak coupling Fermi limit and a new phase ap- 
pearing above U c ~ t for W = 10 — 15. Two numerical evidences of a 
transition are presented. The first is given by the divergence of a char- 
acteristic length of the two dimensional system which allows us to map 
the finite size data onto a single scaling curve. The second is given by 
the existence of a size independent distribution of the first excitation at 
U c . 

In this limit, the use of the parameter r s becomes questionable, since 
it is no longer an appropriate measure of the Coulomb to Fermi energy 
ratio, the one body part of the Hamiltonian being deeply changed by 
the random substrate. The values U c obtained when W is large give 
for the considered filling factors critical values r c s = V r c / '(2£y 'ir(N / 'L' 2 ) 
of the order of the previously obtained ratios rf for a weaker disorder. 
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However, the use of a dimensionless ratio ignoring W when W is large 
can be misleading. 

7.1. DIVERGENCE OF A CHARACTERISTIC 

LENGTH IN TWO DIMENSIONS 
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Figure 21. Right: Characteristic length £z, as a function of L for U < 1 (filled 
symbols): (triangles up), 0.25 (triangles left), 0.5 (triangles down), 0.75 (triangles 
right), 0.85 (circles) and 1 (squares) and U > 1 (empty symbols): 1.15 (squares), 1.25 
(triangles), 1.35 (diamonds) and 1.5 (circles). Left: Ratios £l/L mapped onto the 
scaling curve / as a function of (,/L. The two dimensional scaling length f is given in 
the inset. 



From a finite size study [53], one can obtain scaling laws consistent 
with the divergence of a characteristic length of the two dimensional 
system at a first lower threshold U c ~ t, as in a second order phase 
transition. For defining a characteristic length of the iV-body ground 
state, one considers the change 5pj of the charge density induced by a 
small change 5vi of the random potential Vi located at a distance r = 
\i — j\. To improve the statistical convergence, one computes the change 
Sp(r) = J2j 3pr,jy °f the charge density on the L sites of coordinate j x = 
r yielded by the change vq^ v — ► l.Olvo,^ for the L random potentials 
of coordinate i x = 0. For a Slater determinant made with N occupied 
single particle eigenfunctions (ip a ), fhst order perturbation theory gives: 



Spj = 2Sv, 



N 
a=l /3^a 



^ a {i)^p{i)^ a {j)^fi{j) 



En — E n 



cc cxp 



2r 



(19) 
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the index [3 varying over the one-body spectrum. Therefore, in the 
absence of interaction, the change Sp remains localized on a scale given 
by the one-body localization length (£x, ks L\/2). 

Let us study the dependence of £l on the interaction U for an ensemble 
of 5 x 10 3 clusters, with N = 3, 4, 5 particles in square lattices of size 
L = 24, 28, 31 respectively, corresponding to a very low filling factor 
n e ~ 5 x 10~ 3 . To have Anderson localization inside these sizes we 
considered a large disorder to hopping ratio W/t = 10. Therefore the 
low energy tail of the one body spectrum is made of impurity states 
trapped at some site i of exceptionally low V{. As we are interested in 
studying the effect of Coulomb repulsion on genuine Anderson localized 
states we get rid of the band tail. Typically we ignore the L 2 /2 first 
one-body levels (but results do not change, provided that the Fermi 
level is out of the band tail, ep > — At). From this restricted subset of 
one-body states we built a basis for the iV-body problem, truncated to 
the Nh = 10 3 Slater determinants of lowest energy (convergency tests 
are discussed in Ref. [53]). We note that, for W/t = 10, the one-body 
localization length L\ ~ 4 is smaller than the distance between particles 
d« 15. 

We checked that |<5p(r)| is reasonably fitted by a log-normal distribu- 
tion. Therefore it makes sense to characterize the typical strength of the 
fluctuations by 

fytyp (r) = exp < In | Sp(r) | > (20) 

and extract the length £/, over which the perturbation is effective from 
the exponential decay 

<*Ptyp( r ) « exp(-r/£ L ). (21) 

The size dependence of £l is presented in Fig. 21 left, for increasing 
Coulomb repulsions. One finds the behavior typical of a phase transition: 
£l converges towards a finite value when U < U c (localized phase), 
diverges linearly as a function of L at U = U c ~ 1 (critical point) and 
diverges faster than linearly when U > U c (delocalized phase). This is 
exactly the behavior [54, 55] which characterizes the one-body Anderson 
transition in three dimensions. 

In Fig. 21 right we verify a usual scaling ansatz [54] inspired by the 
theory of second-order phase transitions: 

T = '(?)' <22> 

where we assume that it is possible to map the characteristic length £j, 
at the system size L onto a scaling curve f(L/^), where £ is the scaling 
length characteristic of the infinite two dimensional system. 
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All the data of Fig. 21 left can be mapped onto the universal curve 
/ shown in Fig. 21 right, assuming the scaling length £ given in the 
inset. When U < U c , this length characterizes the localization of the 
effect of a local perturbation of the substrate in the two dimensional 
thermodynamic limit. Our data are consistent with a divergence of £ at 
a threshold U c ~t. 

Very often, additional corrections oc L~ Q to the scaling ansatz occurs 
for small sizes. We point out that our results can be fitted by a simple 
linear law £l = 0.17L for U = U c . This is a further indication that 
the simple ansatz (22) describes scaling for L > 24, without noticeable 
additional L~ a corrections. 

This tells us that the polarized electron system in a highly disordered 
2d substrate becomes correlated when U > U c and that the effect of a 
local perturbation becomes delocalized. 

7.2. SIZE INDEPENDENT DISTRIBUTION 

OF THE FIRST EXCITATION 

In Ref. [56] we used the configuration interaction method, discussed in 
appendix, to study the statistics of the first many-body energy spacing. 
Questions related to the convergence of the method when the Hartree- 
Fock basis is truncated are discussed in Ref. [56]. 

Let us consider the first A-body energy levels Ei, (i = 0,1,2,...) 
for different sizes L, with a large disorder to hopping ratio W/t = 15 
imposed to have Anderson localization and Poissonian spectral statistics 
for the one particle levels at U = when L > 8. This corresponds to 
a strongly disordered limit where the one particle localization length L\ 
is not only smaller than L and the distance d between the particles (as 
in the previous subsection), but becomes also of the order of the lattice 
spacing. We considered A = 4, 9, and 16 particles inside clusters of size 
L = 8, 12, and 16 respectively. This corresponds to a constant low filling 
factor n e = 1/16. We studied an ensemble of 10 4 disorder configurations. 

The first average spacing < AEq > is given in Fig. 22. It exhibits 
a power law decay as L increases, with an exponent [3 given in the 
inset. One finds for the first spacing that (3 linearly decreases from 
d = 2 to 1 when U increases from to 2.5. The next mean spacings 
< AEi >=< Ei + \ — Ei > depend more weakly on U, as shown in Fig. 
22. 

For C/ = 0, the distribution of the first spacing s = AEq/ < AEq > 
becomes closer and closer to the Poisson distribution Pp(s) when L 
increases, as it should be for an Anderson insulator. For a larger U, 
the distribution seems to become close to the Wigner surmise P\y{s) 
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Figure 22. Size dependence of the average gap (first spacing < AEo >oc L ^ u '), for 
W = 15, filling factor n e = 1/16. From bottom to top: [7 = 0, 0.5, 0.8, 1, 1.2, 1.5, 2, 2.5. 
Inset: /3(U) (circles, characterizing < AEo > and squares, characterizing < AEi >, 
with an average over i — 1 — 3). 



characteristic of level repulsion in Random Matrix Theory, as shown 
in Fig. 23 left, for U = 2.5 and L = 16. To study how this P(s) 
goes from Poisson to a Wigner-like distribution when U increases, we 
use the spectral parameter r/o which decreases from 1 to when P(s) 
goes from Poisson to Wigner. In Fig. 23 right, one can see that the 
three curves r]o(U) characterizing the first spacing for L = 8,12,16 
intersect at a critical value U c /t ~ 1.1. Our data suggest that for U < U c 
the distribution tends to Poisson in the thermodynamic limit, while for 
U > U c it tends to a Wigner-like behavior. At the threshold U c , there is a 
size-independent intermediate distribution shown in the inset of Fig. 23 
left, exhibiting level repulsion at small s followed by a exp(— as) decay at 
large s, with a ~ 1.52. Such a size independent distribution is known for 
characterizing a mobility edge in an one body spectrum. This Poisson- 
Wigner transition characterizes only the first spacing, the distributions 
of the next spacings being quite different. The inset of Fig. 23 right 
does not show any intersection for the parameter rj calculated for the 
second spacing. 
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Figure 23. Left: gap distribution P(s) for U — (empty circles) and U — 2.5 
(filled circles) when N — 16, L — 16, W = 15, compared to Pp{s) and Pw{s). Inset: 
size invariant P{s) at U c & 1.1 ; L — 8 (circles), 12 (squares), and 16 (diamonds). 
Right: parameter r]o{U) corresponding to the first spacing ASo at L = 8 (circles), 
12 (squares), and 16 (diamonds). Inset: r]o(U) for the second spacing AEi. 



The observed transition, and the difference between the first spacing 
and the following ones is mainly an effect of the H-F mean field. For the 
first spacings, the curves r]o calculated with the HF data are qualitatively 
similar. At the mean field level the low energy excitations are particle- 
hole excitations starting from the ground state. The energy spacing 
between the first and the second excited states is given by the difference 
of two particle-hole excitations and a Poisson distribution follows if the 
low energy particle-hole excitations are uncorrelated. 

We note that the energy of an electron-hole pair is given by €j — e« — 
U/rij and the classical argument for the existence of a gap in the single 
particle density of states does not apply for the many-body spectrum[57]. 
Therefore the observed opening of a gap for the first energy excitation is a 
remarkable phenomenon beyond the predictions of the classical Coulomb 
gap model. 

7.3. CHANGE OF THE INVERSE 

COMPRESSIBILITY 

In Ref. [56] we studied also the inverse compressibility 

A 2 (N) = E (N + 1) - 2E (N) + E (N - 1), (23) 

i.e. the discretized second derivative of the ground state energy Eq with 
respect to the number N of particles. 
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We consider N = 4, 9, 16 particles on square lattices of size L = 6, 9, 12 
respectively, corresponding to a constant filling factor n e = 1/9. Here we 
focus on the strongly localized regime with disorder strength W = 15. 

The following inverse compressibility data are obtained with the con- 
figuration interaction method. We have checked that the residual inter- 
action does not change qualitatively Hartree-Fock results. Indeed HF 
approximation gives in this weak coupling regime a good estimate of the 
ground state energy (see Ref. [56]) and the inverse compressibility is a 
physical observable which only depends on the ground state energies at 
different number of particles. Neither a precise knowledge of the ground 
state wavefuction (as in the calculation of persistent currents) nor ex- 
cited states energies (as in studies of spectral statistics) are required. 
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Figure 24- Left: Inset: size dependence of the average inverse compressibility < 
A2 >, for W — 15, filling factor n e = 1/9, U = (circles), 0.5 (squares), 1 (diamonds), 
2 (triangles up), and 4 (triangles down). Straight lines are power law fits < A2(U) >oc 
L~ a ^ '. Main figure: exponent a (U). Right: distribution of the normalized inverse 
compressibilities for N — 16, L — 12, W — 15, U — (empty circles) and U — 3 (filled 
circles), fitted by a Gaussian of standard deviation a — 0.30 (dotted line). Dashed 
line gives Poisson distribution. Disorder average is over 10 3 configurations. 



Fig. 24 left shows the L-dependence of the average inverse compress- 
ibility, which is well fitted by the power law < A2(U) >oc L~ a ( u >, with 
a(U) going from 2 to 1 when U goes from to 3 approximately. The 
value a = 2 is expected without interaction (< A2 >=< A >oc 1/L 2 , 
with < A > single particle mean level spacing). The exponent a = 1, 
in this strongly localized regime, can be related to the Coulomb gap in 
the single particle density of states [57]: According to Koopmans' the- 
orem (see Refs. [58, 59, 60] for a thorough discussion about the limits 
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of validity of the Koopmans' theorem for disordered quantum dots), one 
assumes that all the other charges are not reorganized by the addition 
of an extra charge, 

A 2 « ejv+i - e N oc — (24) 

due to the Coulomb gap, where e k is the energy of the k-th. HF orbital 
at a fixed number N of particles. 

The distribution of inverse compressibilities, for L = 12 and W = 15 
is shown in Fig. 24 right. At U = 0, due to Anderson localization, A 2 
distributions are close to the Poisson distribution (deviations from the 
Poisson distribution at small s values are due to the finite system size). 
On the contrary, the distribution at U = 3 shows a Gaussian shape. 
This can be understood within the Koopmans' theorem, since the HF 
energies are given by 



N 

e k = {ipk\Hi\^ k ) + J2 (Q a a k k - Q k £) , (25) 

a=l 

with H\ one-body part of the spinless Hamiltonian and Q^g interaction 
matrix elements given by Eq. (C.3) in appendix C. Due to the small 
correlations of eigenfunctions in a random potential, one can reason- 
ably invoke the central limit theorem (see Ref. [60] for a more detailed 
discussion). We point out that a Gaussian-like distribution of inverse 
compressibilities has been observed in quantum dots experiments in the 
Coulomb blockade regime [61, 62, 63]. This implies that inverse com- 
pressibility fluctuations are dominated by interaction effects instead of 
single particle fluctuations. 

We note that the Gaussian-like shape of the inverse compressibility 
distribution becomes more asymmetric when the disorder strength is 
increased. This asymmetry is due to the Coulomb gap and has been 
discussed in the classical limit [t = 0) in Ref. [64]. 

A simple model for the addition spectra of quantum dots has been 
recently proposed by B. Shapiro [65]. This model explains the change of 
the inverse compressibility distribution when r s increases assuming that 
the system can be decomposed into a stable N-electron Wigner lattice 
plus a lattice of interstitial sites where the added particle can move and 
for which a usual single particle description is assumed. This suggests 
possible relations between the problem of adding an extra particle to a 
Wigner solid and its melting when r s decreases. 
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8. CONCLUSION 

We have numerically revisited the quantum-classical crossover from 
the weak coupling Fermi limit towards the strong coupling Wigner limit, 
using small lattice models. In the continuous limit, the assumed picture 
[17] is relatively simple: a liquid-solid first order transition for r s ~ 37. 
For a clean lattice model at a filling factor 1/9, our results raise the 
question of the possible existence of an intermediate liquid-solid phase 
as first proposed by Andreev and Lifshitz. This may also raise questions 
about the differences between the continuous limit and lattice models. 
When the spin degrees of freedom are included, our lattice model is 
a Hubbard model with additional long range repulsions far from half 
filling. The physics of such models is not very well known, but at least 
suspected to be very rich and complex, as discussed in many works since 
the discovery of high T c superconductors. The role of the disorder is 
another source of complexity. We have divided this chapter in two main 
parts, the first where the disorder is weak, the second where the disorder 
is strong and makes the quantum kinetic effects less relevant. This raises 
also the question of the difference between the quantum Wigner glass 
and the classical Coulomb glass. For large disorder, a finite size scaling 
analysis leads us to conclude that the weak coupling Fermi phase is 
delimited by a second order quantum phase transition, when r s increases. 

Eventually, our first motivation was the experimental discovery of a 
possible new metal for intermediate values of r s . Our numerical studies 
have not directly addressed the transport properties for intermediate r s . 
One cannot claim that the observed intermediate regime gives a new 
metal in the thermodynamic limit. Nevertheless, if transport is mainly 
due to the presence of delocalized defects in a floppy and pinned elec- 
tron solid for intermediate r s , it might be interesting to extend the FLT 
theory a la Landau proposed in Ref. [24] for a clean system to the case 
where one has also elastic scattering by impurities. It will tell us if a 
modified FLT adapted to a very special Fermi system with a non con- 
ventional Fermi surface, still gives rise to usual quantum interferences 
leading to weak localization and eventually to Anderson localization for 
the conjectured gapless excitations of a floppy solid. This could help 
to determine whether one has a true new metal for intermediate r s and 
weak disorder, or only an "apparent" metallic behavior which should 
disappear at the true T — ► limit. This might help to explain why in 
measurements down to 5mK in a GaAs hetrostructure [36], the weak 
localization correction remains less than a few percent of the value pre- 
dicted for a standard disordered 2d Fermi liquid. 
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Appendix: Hartree-Fock approximation 

This is the usual mean field approximation when one has electron-electron inter- 
actions, which we shortly review for the case of spinless fermions. The HF ground 
state (GS) is the Slater Determinant (SD) which minimizes the GS energy expec- 
tation value. In the HF approximation, one reduces the two-body part of the total 
Hamiltonian to an effective single particle Hamiltonian [66, 67, 68] 



u( 52 -r n ^ - 2 t-A c M°*))> (a.i) 



where (...) stands for the expectation value with respect to the HF ground state, 
which has to be determined self-consistently. The first (Hartree) term describes the 
interaction of any electron with the charge distribution set up by all the other elec- 
trons, while the latter (Fock) term is a nonlocal exchange potential. The Hartree 
term comes out as an extra on-site disorder potential, while the Fock term introduces 
extra hopping amplitudes. 

The main advantage of the Hartree-Fock approximation is that it reproduces the 
single-particle density of states g(E) [68], particularly the Coulomb gap near the 
Fermi energy Ep: in the two-dimensional case, g(E) oc \E — Ef\ [57]. The physical 
argument underlying this relation is that an empty site j and an occupied site i 
with a difference in energy Cj — a smaller than 8 must be at a distance larger than 
U/S as the change AEij = ej — e» — U/rij of the system energy must be positive 
when we consider an excitation starting from the ground state. However, the HF 
approximation gives a Coulomb gap also in the delocalized regime at small disorder 
(single particle localization length larger than the system size), where complicated 
many-body effects beyond HF approximation should screen Coulomb interaction. 

Appendix: Stoner instability 

Ferromagnetic instabilities come from the interplay between Coulomb repulsion 
and the Pauli principle. In the Pauli picture, electrons populate the orbital states of 
a system, such as a quantum dot or a metallic grain, in a sequence of spin up - spin 
down electrons. The resulting minimum spin state minimizes the kinetic energy: it 
costs energy to flip a spin since it must be promoted to a higher energy level. On the 
other hand, the maximum spin state requires a maximally antisymmetric coordinate 
wavefunction, thus reducing the effect of Coulomb repulsion. This is at the basis of 
Hund's rule for atoms: electrons occupy orbitals in an open shell so as to maximize 
their total spin. 

Due to the locality of the Pauli principle, ferromagnetic instabilities can be studied 
within the Hubbard Hamiltonian, 

H — -t 22 C L C J> + 2_ / Vim ' 7 + u 2j n iT n a- (B-l) 



In the HF approximation the interaction part of the Hubbard Hamiltonian (B.l) 
is reduced to 



U) n iCT (ni_ CT ). (B.2) 
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The HF energies are given by 

e aa ^ e° a + uJ2\^(i)\ 2 \Mi)f (np-,), (B.3) 

with e% one-body eigenenergies and np- a = d g _ a dp- a occupation numbers for the 
HF orbitals. The total ground state HF energy reads 

E = ^e°(n aCT ) + Uj^ \ip a {i)\ 2 \Mi)\ 2 (n^){n(3i). (B.4) 

In order to study the stability of the nonmagnetic solution, 

; * ea ° < eF (b.5) 

it e aa > e F , 

with t F Fermi level, we take a layer (of thickness 8e) of electrons with spin down 
below the Fermi level to put them in states with spin up [69]. This changes the spin 
state of 8n s = p(e F )Se electrons per unit volume, where p(t F ) is the density of states 
per volume per spin. The change of the one-body energy density is given by 

8w = p(e F )(8e) 2 . (B.6) 

In the clean case, the change of the interaction energy (per volume) is 

*»int = U (y +/»*) (y -P^) - ^ = -U P 2 (e F )(8e) 2 . (B.7) 

Therefore the nonmagnetic state becomes unstable when <5wo + 8w- m t < 0, that is 

Up(e F ) > 1, (B.8) 

which gives the Stoner criterion for ferromagnetic instability. 

In the diffusive regime, the effective interaction strength is enhanced by the pres- 
ence of disorder, leading to ferromagnetic instabilities already below the Stoner thresh- 
old, as pointed out in Ref. [70]. 

We also note that Stoner criterion predicts a ferromagnetic ground state for the 
Hubbard model even at finite temperatures in one and in two dimensions, thus vio- 
lating the Mermin and Wagner's theorem. Therefore, when the Stoner criterion (B.8) 
is satisfied, we can only conclude that the nonmagnetic ground state is unstable. 

The problem of a possible magnetization of the ground state is not only discussed 
in a dilute 2DEG, but is also central in the studies of mesoscopic quantum dots since 
their Ohmic resistances are measured as a function of a gate voltage in the Coulomb 
blockade regime. The possibility of a spontaneous magnetization S of their ground 
state due to electron-electron interactions has been proposed to explain the observed 
conductance peak spacing distributions. 

Appendix: Configuration interaction method 

Even though the approximations involved in the HF method are uncontrolled, the 
mean field HF results can be improved using a numerical method [71, 72, 56] familiar 
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in quantum chemistry as the configuration interaction method (CIM) [73]. Once a 
complete orthonormal basis of HF orbitals has been calculated, 

Hhf(\iPi), ■■-, \lj)N))\lpa} = e a \^p a ), (CI) 

with a = 1,2, ... ,L 2 , it is possible to build up a Slater determinants' basis for the 
many-body problem which can be truncated to the Nh first Slater determinants, 
ordered by increasing energies. The two-body Hamiltonian can be written as 



flint = \ Y. Ql S A4d s d,„ (C.2) 



2 



with 



QlS = U y^ '>Pc{i)'tpl3(j)lp'y{i)lp6U) ,Q g\ 

and djj, = ^2-ipa(j)Cj. One gets the residual interaction subtracting Eq. A.l from 
Eq. C.2. 
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